Generation of arbitrary Dicke states in spinor Bose-Einstein condensates 
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We demonstrate that the combination of two-body collisions and applied Rabi pulses makes it 
possible to prepare arbitrary Dicke (spin) states as well as maximally entangled states by appropriate 
sequencing of external fields. 
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It is now widely recognized that the use of squeezed 
atomic states has the potential for substantial improve- 
ment in the sensitivity of atom interferometers. This 
realization has led to much theoretical and experimental 
work on schemes to realize atomic spin squeezing. Build- 
ing upon the seminal work on photon number squeezing 
by Kitagawa and Yamamoto 0] , Kitagawa and Ueda pro- 
posed and clarified the basic issues concerning the defi- 
nition and preparation of spin squeezed states (SSS's) 
0,|J. This was followed by studies of the production 
of squeezed atomic states [|]|| , the transfer of squeezing 
from incident squeezed light to an assembly of cold atoms 
and the measurement of a collection of squeezed 
atoms |^-[l^|. So far, though, this work seems limited to 
the quantum state control of just a few atoms or ions. 

At the same time, Bose-Einstein condensates have be- 
come almost routinely available in the laboratory. They 
provide us with large coherent ensembles of ultracold 
atoms which seem ideal to perform quantum state manip- 
ulation experiments from an exquisitely well-controlled 
initial state. In particular, recent experiments on two- 
component Bose-Einstein condensates (BECs) in 87 Rb 
p3| have led to considerable work on the dynamics of 
the relative phase and number fluctuations of Bose con- 
densates. There have been vigorous efforts on the issue 
of temporal phase coherence between two coupled Bose 
condensates. A number of workers have thought of this 
system as a prototypical Josephson junction between two 
weakly-coupled superfluid systems. 

While much of the theoretical work so far amounts to 
what are essentially semiclassical (mean-field) analysis 
fTi] |l7t and studies involving small quantum corrections 
going beyond the mean-field Gross-Pitaevskii equation 
jl4],[l8-22|, some publications have gone past this frame- 
work. For instance, Villain et al. |18| , Steel and Collett 
p2fl , and Javanainen and Ivanov |20j ] have studied dy- 
namical aspects of the quantum state of a coupled BEC. 
There has also been intensive work on spin manipula- 
tion in the context of spinor BECs, motivated by the 
experimental work on 23 Na by the group at MIT p3,pi|. 
Theoretical analyses, notably that of Law et al. have 



addressed the issue of collective spin properties in spinor 
BECs |2f| . Recently, S0rensen et al. studied spin squeez- 
ing in a spinor condensate due to the internal nonlin- 
ear atom-atom interaction p6[ . In this note, we extend 
these results by describing a quantum control technique 
which allows to prepare arbitrary Dicke spin states, as 
well as maximally entangled states by means of properly 
sequencing of external coupling fields. 

To set the stage for the discussion and relate it to the 
physics of Bose-Einstein condensates, we first recall some 
basic properties of coherent spin states (CSS's). [^9| Con- 
sider a two-component condensate with iVj atoms in com- 
ponent i, and N1+N2 — N. A general state of the system 
can be written as a superposition of number difference 
states 



I*) = 



N 

£ 

^2=0 



cjv a |JVi,iV 2 



(1) 



Alternatively, we may introduce the effective angular mo- 
mentum quantum number j = (N± + A^)/2 = A^/2 and 
the z-projection quantum number m = (N2 — N{)/2 and 
reexpress the state |^) in terms of angular momentum 
states \j, to). In particular, we have the correspondence 



\j,-j) = \N 1 = N,N 2 = 0) 



(2) 



a state which can be thought of as the "ground state" of 
coherent spin states. Our goal is to construct arbitrary 
states of the system from this ground state \j, —j). 

In a landmark paper, Kitagawa and Yamamoto 
showed that photons could be number-squeezed by mak- 
ing use of a Kerr nonlinear medium. Their procedure 
consisted of splitting a coherent beam, passing one of its 
components through a Kerr nonlinear medium, and inter- 
fering the resulting light with the other coherent beam in 
a Mach-Zchnder interferometer. The enhanced squeezing 
achieved in that case results from the fact that it involves 
a quartic four-wave process, while ordinary squeezing is 
a quadratic process. Following this work, Kitagawa and 
Ueda showed that it was possible to likewise pro- 
duce squeezed spin states (SSS's) by making use of spin 
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Hamiltonians quadratic in the spin operators. Our pro- 
posed scheme relies similarly on the existence of such a 
term in the Hamiltonian describing the dynamics of a 
coupled two-component condensate. 

We consider for simplicity the so-called two-mode 
model describing the coupled two-component condensate. 
This approximation consists in neglecting all modes ex- 
cept the condensate modes. Its second quantized Hamil- 
tonian is 



H = Eiblh + E 2 blk 
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where bi is the annihilation operator for component i, 

Ei = J d 3 r( f>*(r) [p 2 /(2m) + ^(r)] &(r) (4a) 

is the single-particle energy of mode i, V^(r) being the 
trapping potential, and 
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describes the two-body collisions in the condensate in the 
s-wave scattering approximation. Here is the scatter- 
ing length for a two-body collision between an atom of 
« th component and that of j th component and <pi (r) rep- 
resents the condensate wave function for mode i. Finally 



and the real Rabi couplings fl x — Re(il),ft y = Im(fl). 
Note that some control of the nonlinear parameter k can 
be achieved through the proper engineering of the trap- 
ping potential, and hence of the condensate wave func- 
tions 4>i(v). Another way to adjust k is through tuning 
the scattering lengths via Feshbach resonances, although 
this is not always easy to do in practice. However, the 
purpose of the present work is to show that spin squeez- 
ing can be controlled and manipulated with the external 
coupling fields instead of the internal nonlinearity of the 
condensate. In deriving (^), we have assumed that 

E x - E 2 + (N - l/2)(un - u 22 ) = 

a condition that can always be achieved by shifting the 
energy levels of the condensate components. If this con- 
dition is not fulfilled the Hamiltonian ([?]) contains an 
additional term proportional to J z . 

Importantly for the following discussion, we observe 
that the independent temporal control of the Rabi pulses 
characterized by £l x and Q y is a well-established experi- 
mental tool. Depending on the particular system of inter- 
est, these pulses can be in the form of laser light, radio- 
frequency microwaves or magnetic fields. 

With the Hamiltonian (Q) at hand, we now demon- 
strate how an appropriate choice of the external coupling 
fields allows one to generate arbitrary Dicke states from 
the ground state |j, — j). We first observe that a general 
CSS 1 6*, ip) can be created by rotating \ j, —j) by the angle 
9 about the axis n = (siny;, — cos</?, 0), 



n = il d 3 r0i(r)02(r) 



(5) 



describes the strength of the linear coupling between 
components. The procedure to obtain the Hamilto- 
nian (^) and the discussion of the validity of the two- 
mode appro ximation can be found for instance in Refs. 
jl4|,[l8 22 27[. We note that this same Hamiltonian has 
been previously considered by Scott et al. [ p8[ in quan- 
tum studies of self-localization. 

The analysis of Eq. ([|) is greatly simplified by the in- 
troduction of the angular momentum operators 



\6, <p) = exp[-i6>( J x sin ip - J y cos <p)] \j, -j) 



(9) 
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Jv = \{b\b 2 -blb x ), 
J z = ^(blb 2 -b\h), 



(6a) 
(6b) 
(6c) 



in terms of which Hamiltonian (|^) may be rewritten as 

H = — £l X J X Qyjy, (7) 

where we have introduced the effective nonlinear coupling 



K = 2^ Uu + U22 ' ~ Ul2 ' 



(8) 



Starting from state switching on the coupling 

pulses and properly adjusting the strength of fl x and Q y , 
we can prepare any CSS. For the CSS \0,<p), the expec- 
tation value and the variance of operator J z are 

(Jz) = ~j cos 6, 
<AJ 2 ) = |sin 2 0. 

(Note that neither (J z ) nor (A J 2 ) will change after the 
coupling fields are turned off.) Following the usual con- 
vention we say that a j-spin state is squeezed along J z if 
the state has the same ( J z ) but reduced (A j 2 ) compared 
to the CSS \0, f). Or we can define a parameter |35| 
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A state is squeezed along j z if £ z < 1. 

Procedures to achieve spin squeezing in the equato- 
rial plane (J z ) = was previously proposed by Kitagawa 
and Ueda and more recently by Law et al. fl30| . 

We briefly recall them as a preparation for our general 
scheme. In the first scheme, a pulse is applied on the 
ground state \j, —j) along the y-direction at t = 0. The 
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duration of the pulse is assumed to be short enough that 
the effects of the nonlinear interaction are negligible dur- 
ing this time. This pulse aligns the spin vector along the 
negative rc-axis, with the resultant state, \J,—J) X . Af- 
ter the pulse, the system evolves under the Hamiltonian 
H sp in — kJ 2 until a second short pulse is applied along 
x at t = t%. At the end of this pulse, we have 



(Jz 
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<A^) = £ 1 
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A- VA 2 + B 2 cos2(a + S) }, (lib) 



where a = — j Q, x dt/h is the strength of the second pulse 
and 8 = \ tan _1 (S/.4) with 

A = l-cos 2(j - 1) (2Ki 1 ), 

B = 4sin(2a) cos 20-1 ^^) sin(/rti). 

Optimal squeezing in J z is achieved for a = —5 and 
\ K \t! = (3/8y/ 6 j- 2 / 3 , which leads to (AJ 2 ) mm « 

l/i\l/3 

2 V3/ . — . 

As pointed out by Law et al. Ref. 1301, an alternative 
approach can be realized by switching on a coupling field 
along the x-direction immediately after the first ir/2 pulse 
along the y-direction, so that the system evolves under 
the Hamiltonian H = nJ 2 — fl x J x . During the evolution 
following the control pulses, ( J z ) remains equal to while 
(A J 2 ) can be significantly reduced at appropriate time. 
Fig. pi shows the time evolution of the parameter £ z for 
both pulse sequences. 



states with arbitrary prescribed value of (J 2 ) = too and 
arbitrary Dicke states, we propose to use instead the fol- 
lowing method: First, a coupling pulse of appropriate 
strength is applied to the ground state \j,—j), creating 
an initial CSS with (J z ) = m . Assume without loss of 
generality that that pulse is along the x-axis. Following 
this pulse, we then apply a coupling field along the nega- 
tive x-axis. During the subsequent time evolution, both 
( J z ) and (AJ 2 ) start to oscillate. But at certain times 
when (J z ) returns back to its initial value mo, (AJ 2 ) 
comes close to a local minimum which is less than the 
initial variance, i.e., £ z < 1. 
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FIG. 1. Time evolution o f th e parameter £ 2 for j = 500. 
Curve 1 corresponds to Eq. (Ill ) with a — —8; the other two 
curves are obtained using the scheme of Ref. [30] . For curve 
2, Q x /k, = 20; for curve 3, Q x /n — 10. 



The essential point of both schemes is as already men- 
tioned that ( J z ) remains equal to zero after the first 7r/2 
pulse. This, however, needs not to be the case: If the first 
pulse creates a CSS with (J z ) = m ^ 0, then following 
the above procedures one might obtain a state squeezed 
along J z but with ( J z ) ^ to . In order to create squeezed 



FIG. 2. (a) Time evolution of (j z ) (left panel) and £ z 
(right panel). The initial CSS is \0 = Tv/3,ip = 0) with 
j = 500. After the initial CSS is created, the system evolves 
under the Hamiltonian H — kJ z — Q X J X with Q, x /k = —30. 
At times indicated by the arrows, (J z ) reaches the initial 
value (-250) while £ z is reduced, (b) Probability distribu- 
tion of the azimuthal number m. Dashed curve: the ini- 
tial CSS \8 — tt/3, tp = 0); solid curve: the spin squeezed 
state at Kt = 0.0344. (c) Quasiprobability distribution on the 
Bloch sphere for the initial CSS and the squeezed state at 
Kt = 0.0344. 



Provided that the coupling field is turned off at these 
precise times, a state squeezed along J z with (J z ) = too 
will be created. For large enough squeezing, the state 
thus prepared can be regarded as an approximation of 
the Dicke state \j, too) which is the spin analog of the 
Fock state. A typical time evolution of (J z ) and £ z is 
shown in Fig. |](a) where the initial CSS has j = 500 
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and (j z ) — —250. In this example, the optimal squeez- 
ing occurs at nt = 0.0344 where the variance (A Jf) is 
reduced by a factor of more than 20. Fig. ||(b) shows the 
probability distribution against the azimuthal number m 
at the optimal squeezing, while in Fig. ||(c) we plot the 
quasiprobability distribution (QPD) on a Bloch sphere. 

Squeezing and entanglement are oftentimes closely re- 
lated. In particular, it is well known that the quadratic 
Hamiltonian H SV m = kJ% also produces entanglement 
in a collective spin system. In particular, starting from 
the ground state \j,—j) x , H sp i n generates at time t* = 
Tin /(2k) the state 
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|j,-j) x + (-l)i e W4| i)+i ) a . (12) 



The state is called maximally entangled because if one 
individual spin is found to be aligned along the nega- 
tive x-axis or the positive x-axis, so are all other spins. 
However, the degree of entanglement degrades after t*, 
a consequence of the fact that the state ( |l2|) is not an 
eigenstate of the Hamiltonian H sp i n . To preserve the 
entanglement in our system, a 7r/2 pulse along the y- 
direction can be applied at t = t*. This converts the 
state \j,±j) x into \j,±j). The resultant state, still max- 
imally entangled, it now an eigenstate of H sp i„, hence 
it remains maximally entangled during the subsequent 
time evolution. This state is also a Schrodinger cat state 
since it is a coherent superposition of two macroscopi- 
cally distinct states — one state has all the population 
in component 1 and the other in component 2. 

Maximally entangled states, in particular those of mas- 
sive particles instead of fast-escaping photons, are of 
great importance in fundamental physics as well as in 
applications in quantum information and quantum mea- 
surement. A great deal of effort has been directed to- 
ward the creation of such states. Two- |3^j3^], three- 
pi^ and four-particle |34|] entanglement have been suc- 
cessfully demonstrated experimentally in trapped ions, 
Rydberg atoms, and cavity QED. However, a further in- 
crease of the number of entangled particles in these sys- 
tems are expected to be a severe experimental challenge. 
The two-component, atomic BEC with its built-in intrin- 
sic nonlinearity appears to be a promising candidate to 
generate entanglement on a macroscopic level. 

To conclude, we have shown that an arbitrary collec- 
tive spin squeezing and entanglement of a two-component 
spinor condensate can be readily controlled by the cou- 
pling fields between the two components. Squeezed or en- 
tangled spin states will find applications in high-precision 
spectroscopy, atomic interferometry and quantum infor- 
mation, and spinor condensates are attractive candidates 
to create such states. In practice, the scheme presented 
in this work can be realized in the two-component 87 Rb 
condensate 13 . Another possibility is offered by the 
F = 1 spinor condensate of 23 Na |23 24] which can be 



reduced to an effective spin-(l/2) system by shifting the 
energy level of the mp = state, e.g., with an ac Stark 
shift. 
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